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Abstract 

A family of quantum anharmonic oscillators is studied in any finite spatial dimen- 
sion in the scheme of first quantization and the investigation of their eigenenergies is 
presented. The statistical properties of the calculated eigenenergies are compared with 
the theoretical predictions inferred from the Random Matrix theory. Conclusions are 
derived. 

1 Motivation 

The quantum harmonic oscillator proved to be a fructuous model of many physical systems: 
quantum electromagnetic field or systems of atoms (ions, nuclei) in ideal crystals interacting 
via harmonic attractive force, etc. In the former case the excitation particles or quanta of 
the electromagnetic field are called photons [1, 2] whereas in the latter case the elementary 
excitation particles of vibrations of crystal lattice or quanta of the sound field are named 
phonons [3, 4]. Both of these quantum fields are bosonic ones [5, 6]. Also in the case of 
the interaction of the quantum electromagnetic field with the matter field via dipolar elec- 
trostatic interaction the quantum harmonic oscillator is hugely investigated. The harmonic 
potential energy is only an approximation for the real anharmonic potential energy of mu- 
tual interaction between atoms (ions, nuclei) in real crystals. Therefore the motivation of 
the present work is a more realistic description of quantum anharmonical systems. 



2 Quantum harmonic oscillator in D = 1 spatial dimen- 
sion 

Firstly: Our study commences to concentrate on the case of simple quantum harmonic 
oscillator in D = 1 spatial dimension. The dimensionless Cartesian coordinate is denoted 



1 



by Xj and its conjugated dimensionless linear momentum is pj,j > 1. Let us consider the 
Hilbert space: 

Vi = L2(R,C,dx,), (1) 

of the complex-valued wave functions that are modulus square integrablc on the set R 
of the real numbers, and C is the set of the complex numbers. The Hilbert space Vi is 
a separable space and its orthonormal basis Bi is a set of Hermitc's functions (Fock's 

functions, eigenfunctions of the dimensionless Hamiltonian operator Ti^j^ of the quantum 
harmonic oscillator in D — 1 spatial dimension): 

^r,{xj) = Nr,Hr,{xj)eM-lx'),Nr, = [V^rj\2'TKrj e N, (2) 
where N is a set of natural numbers including zero, whereas: 

Hrj{xj) = {-ly^ exp(a;|)^^ exp(-a;|), (3) 

is Tj-th Hermite's polynomial [7]. We assume from now that all the investigated operators 
are dimensionless (nondimensional) . The dimensionless quantum operator Xj of the jth 
component of the position vector (radius vector) operator x acts on the Hermite's basis ^^^^ 
function as follows [7]: 

= \/f + (4) 

whereas the quantum operator pj of the j'th component of linear momentum vector operator 
p is given in the basis by the following formula: 

Matrix elements {xj)i-^rj and {pj)i.^r. of these operators equal correspondingly: 



/oo 
%{Xj)Xj^!r^{Xj)dXj = 
-oo 



where 
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is discrete Kronecker's delta (it is not continuous Dirac's delta distribution 5-d) . The quantum 
operator x| of the square of the operator Xj and quantum operator p| of the square of the 
operator pj have the following representation in the basis Bi. 

xJ*,, = |[v/^7^^*,._2 + {2rj + 1)*,, + ^/¥-+T,/¥-+l^r,+2l (9) 
Pprj = |[-^A7^^*.,-2 + (2r,- + 1)*,. - V^7TTv/^7T^*.,+2], (10) 

and their matrix elements {xj)i.^rj, ip'j)ij,rj read: 

ip')h,r, = ^[-^r,-1^6i^,r,-2 + (2r, + 1)5,^,,^. - ^rj + l^r, + 26i^,r,+2]. (12) 

The dimensionless (nondimensional) quantum Hamiltonian operator 'H(j) of the quantum 
harmonic operator in one spatial dimension is defined as follows: 

The basis function ^^^^ is its eigenfunction: 

7Y(,)*,, = e,.*,,, (14) 
therefore its matrix element ('^(j))/j,rj is equal: 

i^{3))lj,rj ^^rjSl.^rj, (15) 

where 

e,, = 2rj + 1, (16) 

is the Tj-th eigenenergy of 'H(j). The eigenenergies are simply all odd natural numbers, and 
the quantum Hamiltonian is diagonal operator, and its matrix representation is diagonal 
CX3 X CX3 matrix. Note, that if one introduces the notion of nearest neighbour energy spacing 
(NNS) Srj between two adjacent eigenenergies: 

Srj — ^Tj+l ~ ^rj, (17) 

then for the quantum harmonic oscillator it holds: 

— 2 — const, (18) 

so the eigenenergies are equidistant. Let us consider first N consecutive energy levels e^^ , = 
0, ■ • • , [N—l), N > 2. One spans over the set of the eigenvectors ) '"j = 0, • ■ ■ , (A^— 1), N > 
2, a Hilbert space Vi^n that is a subspace of the Hilbert space Vi. The truncated Hilbcrt space 
Vi,7v is isomorphic to the A^th Cartesian product of complex spaces C: Vi^at = C^. Here 
iV is a complex dimension of truncated Hilbert space Vi^n and of C^. The eigenfunctions of 



iii 



are A^-component constant complex vectors (analogs of A^-spinors), and the operators 
acting on it are N x N deterministic complex-valued matrices. The probability distribution 
Pjv-i of the spacings is discrete one point distribution for any finite value of N, N > 2: 

tending in the thermodynamical limit ^ oo to the singular Dirac's delta distribution 6^: 

Poo(s) = 5d(s-2). (20) 

Secondly, let us perform more difficult task consisting of calculating all the fiip-fiop 
transition amplitudes (hopping amplitudes) from the quantum state Xr] = ^j^^r^ to the 
quantum state {sj > 0): 

/oo 
%{xj)x'/<lfrA^j)dXj. (21) 
-oo 

Physically, the transition amplitude {'ms.)i.^rj is connected with the processes of emissions 
and/or absorptions of Sj phonons, because: 

xy=V2'^(a, + a+)^^ (22) 

where the bosonic single phonon annihilation and creation operators in one spatial 

dimension, respectively, and: 



aj-^r, = v^^r,-i, aj-^r, = sjTj + l^'r.+l- (23) 

One can calculate the lower transition amplitudes manually, e. g., using recurrence relations, 
matrix algebra, etc., but it is tedious (even for 3 < Sj < 6). If one wants to calculate aJi 
the transition amplitudes then he must return to the beautiful XIX century mathematics 
methods and after some reasoning he obtains the exact formula: 

(^Sj)ij,rj 



A,=0 K,=0 



A,!(/,-2A,)! Kj\{Tj-2Kj)\ 



.2f +3-2A,-2«i-i . p^ ^i + ^i + ^j ^^i + ^ )]^ (24) 

where [•] is entier (step) function, F is Euler's gamma function (compare our result Eq. (24) 
with the formulae in [8, 9]). Therefore, the matrix representations of the even power operators 
x'^p in the basis Bi are hermitean (symmetrical real) matrices with nonzero diagonal and 
nonzero Pj subdiagonals (and nonzero pj superdiagonals), where the distance of the nearest 
superdiagonals (or subdiagonals) is 2 (the diagonal is also distant by 2 from the nearest super- 
and sub-diagonal), whereas the odd power operators xj'^^^ in the basis Bi are hermitean 
(symmetrical real) matrices with zero diagonal and pj nonzero subdiagonals (and pj nonzero 
superdiagonals), where the distance of the nearest superdiagonals (or subdiagonals) is 2 (the 
nearest super- and sub-diagonal are also distant by 2). The physical interpretation of the 
superdiagonals (or subdiagonals) is connected with the absorption (or emission) of phonons. 
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3 Quantum anharmonic oscillator in D = 1 spatial di- 
mension 



Thirdly, we are ready to deal with the quantum anharmonic oscillator in D — 1 spatial 
dimension. Its dimensionless Hamiltonian operator 7i! anharm reads: 



(j) ,axiharm. 

Si=0 



^m.a^harm = ^0) + E ^^j^V ' (25) 



where Sj is a degree of the anharmonicity of the oscillator, and the prefactors a^^. are the 
strengths of anharmonicity. The matrix elements of the anharmonic Hamiltonian operator 
are: ^ ^ 

Sj=0 Sj=0 

where the representation of the quantum anharmonic oscillator in the quantum harmonic 

oscillator basis Bi is mathematically correct, because the basis Bi is a complete set, and the 
Hilbert space of the eigenfunctions of the anharmonic oscillator is isomorphic to the Hilbert 
space Vi for the harmonic oscillator, provided that the total potential energy U^^^ ^^^g^j of the 
anharmonic oscillator: 

^/),total(^j) = + ^/),anharm(^i)' (27) 

is bounded from below (there are no scattering eigenstates) , where the anharmonic potential 
energy U^)^^^ reads: 

^(lanha.m(^.) = E <^s,^ ' (28) 
s,=0 

It suffices that the degree of the anharmonicity Sj = 2S'j is an even number and that the 

strength of anharmonicity as^ is strictly positive: > 0, so that tota\i-'^j) ~^ 
\Xj\ — > oo. 

Fourthly, we repeat the "Bohigas conjecture" that the fluctuations of the spectra of the 
quantum systems that correspond to the chaotic systems generally obey the spectra of the 
Gaussian random matrix ensembles. The qiiantum intcgrable systems correspond to the 
classical integrable systems in the semiclassical limit [10, 11]. The probability distributions 
Pp of the nearest neighbour spacing for the Gaussian orthogonal ensemble G0E(2) of 2 x 
2 Gaussian distributed real-vahicd symmetric random matrix variables (/3 = 1), for the 
Gaussian unitary ensemble GUE(2) of 2 x 2 Gaussian distributed complex-valued hermitean 
random matrix variables (/? = 2), for the Gaussian symplectic ensemble GSE(2) of 2 x 2 
Gaussian distributed quaternion- valued selfdual hermitean random matrix variables (/3 = 4), 
and for the Poisson ensemble (PE) of the random diagonal matrices with homogeneously 
distributed eigenvalues on the real axis R are given by the formulae: 

Pp{s) = e{s)Af,s^exp{-Bps^), (29) 
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for the Gaussian ensembles, and 



Pois) = e{s)cxp{-s), (30) 

for the Poisson ensemble, where 9 is Heaviside's unit step function [12, 13, 14, 15, 16, 17, 18, 
19, 20, 21]. The constants: 

_ r/^+H(/3 + 2)/2) _rW + 2)/2) 

^^-^r^+2((^ + i)/2) ^''-r2((/?+i)/2)^ (^^^ 

arc given by the formulae: Ai = n/2, Bi = 7r/4 (GOE), A2 = B2 = i/n (GUE), 

and A4 = 262144/7297r'^, B4 = 64/97r (GSE), respectively. For the Gaussian ensembles 
the energies are characterized by the "level repulsion" of degree /3 near the origin (near the 
vanishing spacing s = 0), and the probability distributions P/s vanish at the origin, and 
the quantum system with "level repulsion" is cast to the class of quantum chaotic systems. 
For the Poisson ensembles the energies are characterized by the "level clustering" near the 
origin, and the probability distributions Pq has maximum at the origin, and the quantum 
system with "level clustering" are treated as the quantum integrable system. After many 
numerical experiments conducted with different quantum anharmonic oscillators (up to the 
sextic quantum anharmonic oscillator Sj = 6) we draw conclusion that majority of them 
behaves like quantum integrable systems, the eigenenergies tend to cluster, the histogram of 
nearest neighbour spacing is closer to the Pq distribution resulting from the Poisson ensemble 



4 Quantum harmonic oscillator in D > 1 spatial dimen- 
sions 

Fifthly: Quantum harmonic oscillator in D spatial dimensions is a solvable analytically 
model. In order to make the deliberations easier we reduce our present interest to the first 
quantization case. Therefore the relevant Hilbert space Vd is isomorphic to a D-dimensional 
tensor (Cartesian) product of the one- dimensional Hilbert spaces Vi: 

D 

Vd = (8)Vi, (32) 

whereas the rth harmonic oscillator's eigenfunction in D dimensions is, neither sym- 
metrized nor antisymmetrized, tensor product of the eigenfunctions in one dimension: 

D 

*r(x) = l[^rjixj),ic= {xi,...,Xd) G R^, r = (n , . . . , t^,) G N^, (33) 
where we used boldface font for the £)-tuples x, and r. It follows that: 

D D 

^'r(x) = NrH^i^) exp(--x2), iVr = n ^r,, H,{^) = [] Hr.{Xj), (34) 

vi 



whereas 

D 



E(^.)'- (35) 



One can also draw a conclusion that the Hilbert space: 

Vd = L2(R^,C,dx),D > 1, (36) 

is composed of the complex-valued wave functions ^ that are modulus square integrable 
on the set R^. The Hilbert space Vd is separable space, and its orthonormal basis Bd 
is a set of Hermite's functions \&r in D dimensions (Fock's functions in D dimensions). 
The dimensionlcss (nondimensional) quantum Hamiltonian operator Tio of the quantum 
harmonic oscillator in D dimensions is a sum of quantum Hamiltonian operators of the 
quantum harmonic oscillators in one dimension: 

no^f: Hu) = E(p? + %') = + (37) 

The quantum Hamiltonian operator Hd is diagonal in Fock's basis of its eigenfunctions ^^r^ 

Hz)^r = er*r, (38) 

and its matrix element {Tio^r is equal: 

{nD)x,r = er5l,r, (39) 

whereas ^ ^ 

er = Ee'-. =E(2r. + l)' (40) 

J=l 3=1 

is the rth eigenenergy of TY^i and where 

Kr=f{h,r,. (41) 

i=i 

is discrete Kronecker's delta in D dimensions (it is not continuous Dirac's delta in D di- 
mensions). The eigenenergies er are simply the sums of all odd natural numbers, and the 
quantum Hamiltonian TY^i is (direct) sum of diagonal operators in one dimension, and its ma- 
trix representation is direct sum of diagonal oo x oo matrices (it is poly- index matrix) . There 
is no unique straightforward analog of the nearest neighbour spacing in D dimensions, 
for D>2. 

Sixthly, let us perform very difficult task consisting of calculating all the flip-flop transi- 
tion amplitudes (hopping amplitudes) from the quantum state = x^^r to the quantum 
state *i (s > 0): 

P roo 

(ms)l,r = (X^)l,r = (*l|i^*r)v. = \{ / %{x ^)x'^ ^> r,{x ^)dx ^ . (42) 
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The flip-flop transition amplitude (ms)i,r is connected with the processes of emissions and/or 
absorptions of s phonons in D spatial dimensions, because: 

= Ili^?) = fli^^'i^^ + ^tP)] = [V2(a+ a+)]^ (43) 

where a = (ai, ...,&£)), a+ = {af, ...,a~^), are the bosonic multiphonon (D-phonon) annihila- 
tion and creation operators in D spatial dimensions, respectively, and: 

= V^*(ri,...,r-,-l,...,ro)> «/*r = ^r^' + Wo) ■ (44) 

It can be easily proven that: 

(X^)l,r=n(^?)^„..- (45) 

One can calculate the lower transition amplitudes manually, e. g., using recurrence relations, 
matrix algebra, etc., but it is tedious (even for 3 < Sj < 6). The exact formula for all the 
transition amplitudes reads: 

D 

{ms)l,r = l[imsj)lj,r^- (46) 

5 Quantum anharmonic oscillator in D > 1 spatial di- 
mensions 

SeventMy, we are ready to investigate the quantum anharmonic oscillator in D spatial di- 
mensions. Its dimensionless Hamiltonian operator Tif) anharm reads: 

S (Si,...,Sd) D 

^Aanharm = + E = + E hsi,..,sn) ' (47) 

s=0 (si,...,sd)=(0,...,0) i=l 

where S is a D-tuple of degrees of the anharmonicity of the oscillator, and the prefactors 
Og are the strengths of anharmonicity. The matrix elements of the anharmonic Hamiltonian 
operator are: 

s s 

(^D,anharm)l,r = Crt^l.r + E «s(x^)l,r = Crf^l.r + E «s(ms)l,r, (48) 

s=0 s=0 

where the representation of the D-dimensional quantum anharmonic oscillator in the quan- 
tum harmonic oscillator basis Bd is mathematically correct, because the basis Bd is a com- 
plete set, and the Hilbert space of the eigenfunctions of the anharmonic oscillator is isomor- 
phic to the Hilbert space Vd for the harmonic oscillator, provided that the total potential 
energy ^;otai quantum anharmonic oscillator in D dimensions: 

WB,total (X) = + WB,anharm (x) , (49) 
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is bounded from below (there are no scattering eigenstates in D dimensions), where the 
anharmonic potential energy W^^anharm i^" 

s=0 (si,..,sd)=(0,..,0) J=1 

It suffices that the D-tuple of degrees of the anharmonicity of the oscillator S = 2S' is 
composed of even numbers and that the strength of anharmonicity as is strictly positive: 
as > 0, so that ^^^^^^^[(x) oo for |x| oo. 

Eighthly, we repeat that the "Bohigas conjecture" also holds for the quantum oscillators 
in D dimensions. Having conducted many numerical experiments with different quantum 
anharmonic oscillators (up to the sextic {D — 3)-dimensional quantum anharmonic oscilla- 
tors: Sj = 6) we draw conclusion that some of them behave like quantum integrable systems, 
the eigenenergies tend to cluster, the histograms of nearest neighbour spacing are closer to 
the Pq distribution resulting from the Poisson ensemble, whereas other ones look like quan- 
tum chaotic systems, their eigenenergies are subject to repulsion, the histograms of NNS are 
closer to the distributions Pi, P2, -P4, derived from the Gaussian Random Matrix ensembles 
[22]. 
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